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Abstract We discuss the high-density nuclear equation of state within the Brueckner- 
Hartree-Fock approach. Particular attention is paid to the effects of nucleonic 
three-body forces, the presence of hyperons, and the joining with an eventual 
quark matter phase. The resulting properties of neutron stars, in particular the 
mass-radius relation, are determined. It turns out that stars heavier than 1.3 solar 
masses contain necessarily quark matter. 
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1. Brueckner theory 

Over the last two decades the increasing interest for the equation of state 
(EOS) of nuclear matter has stimulated a great deal of theoretical activity. Phe- 
nomenological and microscopic models of the EOS have been developed along 
parallel lines with complementary roles. The former models include nonrel- 
ativistic mean field theory based on Skyrme interactions [1] and relativistic 
mean field theory based on meson-exchange interactions (Walecka model) [2]. 
Both of them fit the parameters of the interaction in order to reproduce the em- 
pirical saturation properties of nuclear matter extracted from the nuclear mass 
table. The latter ones include nonrelativistic Brueckner-Hartree-Fock (BHF) 
theory [3] and its relativistic counterpart, the Dirac-Brueckner (DB) theory [4], 
the nonrelativistic variational approach also corrected by relativistic effects [5], 
and more recently the chiral perturbation theory [6]. In these approaches the 
parameters of the interaction are fixed by the experimental nucleon-nucleon 
and/or nucleon-meson scattering data. 

For states of nuclear matter with high density and high isospin asymmetry 
the experimental constraints on the EOS are rather scarse and indirect. Differ- 
ent approaches lead to different or even contradictory theoretical predictions 
for the nuclear matter properties. The interest for these properties lies, to a 
large extent, in the study of astrophysical objects, i.e., supernovae and neu- 
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tron stars. In particular, the structure of a neutron star is very sensitive to the 
compressibility and the symmetry energy The neutron star mass, measured 
in binary systems, has been proposed as a constraint for the EOS of nuclear 
matter [7]. 

One of the most advanced microscopic approaches to the EOS of nuclear 
matter is the Brueckner theory. In the recent years, it has made a rapid progress 
in several aspects: (i) The convergence of the Brueckner-Bethe-Goldstone 
(BBG) expansion has been firmly established [8, 9]. (ii) Important relativis- 
tic effects have been incorporated by including into the interaction the virtual 
nucleon-antinucleon excitations, and the relationship with the DB approach 
has been numerically clarified [10]. (iii) The addition of microscopic three- 
body forces (TBF) based on nucleon excitations via pion and heavy meson 
exchanges, permitted to improve to a large extent the agreement with the em- 
pirical saturation properties [10-12]. (iv) Finally, the BHF approach has been 
extended in a fully microscopic and self-consistent way to describe nuclear 
matter containing also hyperons [13], opening new fields of applications such 
as hypernuclei [14] and a more realistic modeling of neutron stars [15, 16]. 

In the present article we review these issues and present our results for neu- 
tron star structure based on the resulting EOS of dense hadronic matter. 

Convergence of the hole-line expansion. The nonrelativistic BBG expan- 
sion of the nuclear matter correlation energy E/A can be cast as a power series 
in terms of the number of hole lines contained in the corresponding diagrams, 
which amounts to a density power expansion [3]. The two hole-line truncation 
is named the Brueckner-Hartree-Fock (BHF) approximation. At this order the 
energy D2 is very much affected by the choice of the auxiliary potential, as 
shown in Fig. 1 (solid lines), where the numerical results obtained with the 
gap and the continuous choice are compared for symmetric nuclear matter as 
well as neutron matter. But, as also shown in the same figure [8, 9], adding 
the three-hole line contributions D3, the resulting EOS is almost insensitive to 
the choice of the auxiliary potential, and very close to the result D2 with the 
continous choice. 

In spite of the satisfactory convergence, the saturation density misses the 
empirical value po = 0.17 fm~ 3 extracted from the nuclear mass tables. This 
confirms the belief that the concept of a many nucleon system interacting with 
only a two-body force is not adequate to describe nuclear matter, especially at 
high density. 

Relativistic corrections. Before the possible effects of TBF are examined, 
one should introduce relativistic corrections in the preceding nonrelativistic 
BHF predictions. This is done in the Dirac-Brueckner approach [4], where 
the nucleons, instead of propagating as plane waves, propagate as spinors in 
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Figure 1. Comparison of BHF two hole-line (lines) and three hole-line (markers) results for 
symmetric nuclear matter (left plot) and pure neutron matter (right plot), using continuous and 
gap choice for the single-particle potentials. 



a mean field with a scalar component Us and a vector component Uy, self- 
consistently determined together with the G-matrix. The nucleon self-energy 
can be expanded in terms of the scalar field, 

E(p) = y / (M + Us) 2 +p 2 + U v ^ Po + U v + — Us + ^ 1 U 2 S + ... . (1) 

Po ^Pq 

The second-order term can be interpreted as due to the interaction between two 
nucleons with the virtual excitation of a nucleon-antinucleon pair [17]. This 
interaction is a TBF with the exchange of a scalar (a) meson, as illustrated by 
the diagram (d) of Fig. 2. Actually this diagram represents a class of TBF with 
the exchange of light (tt, p) and heavy (a, u) mesons. There are, however, sev- 
eral other diagrams representing TBF, Fig. 2(a-c), which should be evaluated 
as well in a consistent treatment of TBF. 

2. Three-body forces 

Since long it is well known that two-body forces are not enough to explain 
some nuclear properties, and TBF have to be introduced. Typical examples 
are: the binding energy of light nuclei, the spin dynamics of nucleon-deuteron 
scattering, and the saturation point of nuclear matter. Phenomenological and 
microscopic TBF have been widely used to describe the above mentioned prop- 
erties. 

In the framework of the Brueckner theory a rigorous treatment of TBF 
would require the solution of the Bethe-Faddeev equation, describing the dy- 
namics of three bodies embedded in the nuclear matter. In practice a much sim- 
pler approach is employed, namely the TBF is reduced to an effective, density 
dependent, two-body force by averaging over the third nucleon in the medium, 
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Figure 2. Various diagrams contributing to the microscopic TBF. 



taking account of the nucleon-nucleon correlations by means of the BHF defect 
function gij, 

(l2\V(p)\l'2') = *123 <123|y|l'2 / 3 , > V m . (2) 

33' 

Here ^123 = (/>3(1 — (713) (1 — 523) and fa is the free wave function of the third 
particle. This effective two-body force is added to the bare two-body force and 
recalculated at each step of the iterative procedure. 

Microscopic TBF. The microscopic TBF of Refs. [10-12] is based on meson- 
exchange mechanisms accompanied by the excitation of nucleonic resonances, 
as represented by the diagrams plotted in Fig. 2. Besides the TBF arising from 
the excitation of a NN pair [diagram (d)], already discussed in the preceding 
section, another important class of TBF [diagram (a)] is due to the excitation 
of the isobar A(1232) resonance via the exchange of light (tt, p) mesons, or 
the lowest non-isobar nucleon excitation iV*(1440) excited by heavy meson 
{a and oj) exchanges. Diagrams (b) and (c) are included only for completeness 
and play a minor role [10]. 

The combined effect of these TBF is a remarkable improvement of the satu- 
ration properties of nuclear matter [12]. Compared to the BHF prediction with 
only two-body forces, the saturation energy is shifted from —18 to —15 MeV, 
the saturation density from 0.26 to 0.19 fm~ 3 , and the compression modulus 
from 230 to 210 MeV. The spin and isospin properties with TBF exhibit also 
quite satisfactory behaviour [18]. 
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Phenomenological TBF. A second class of TBF that are widely used in the 
literature, in particular for variational calculations of finite nuclei and nuclear 
matter [5], are the phenomenological Urbana TBF [19]. We remind that the Ur- 

r2w , 



bana IX TBF model contains a two-pion exchange potential V^l supplemented 

rR 
ijk'' 



by a phenomenological repulsive term V- R 



ijk ' ij 

where 



Vijk - Vjjl + V? k , (3) 



eye 



{Xij,Xj k } {Ti-Tj,Tj-Tk} + - [Xij,Xjk] [Ti-Tj,Tj-Tk] 



(4) 



V W = uY, T2 { m ^nj)T 2 {m,r ]k ) . (5) 

eye 

The two-pion exchange operator Xij is given by 

Xij = Y(m n rij)(Ti-crj + T(mT T r ij )S i j , (6) 

where SV, = 3(<Tj -r^^Oj -r^) — <Tj (Tj is the tensor operator and a and r 
are the Pauli spin and isospin operators. Y and T are the Yukawa and tensor 
functions, respectively, associated to the one-pion exchange [19]. 

After reducing this TBF to an effective, density dependent, two-body force 
by the averaging procedure described earlier, the resulting effective two-nucleon 
potential assumes a simple structure, 



Vt n0 (r) 



(vi-aAVFW + SijWVFir) +V H (r), (7) 



tRi 



containing central and tensor two-pion exchange components as well as a cen- 
tral repulsive contribution. For comparison, the averaged microscopic TBF 
[10] involves five different components: 



+ S tJ (r) { Tl . T3 )Vf{r) + V T {r) 



(8) 



In the variational approach the two parameters A and U are determined by 
fitting the triton binding energy together with the saturation density of nuclear 
matter (yielding however too little attraction, E/A « —12 MeV, in the lat- 
ter case [5]). In the BHF calculations they are instead chosen to reproduce the 
empirical saturation density together with the binding energy of nuclear matter. 
The resulting parameter values are A = -0.0293 MeV and U = 0.0048 MeV 
in the variational Urbana IX model, whereas for the optimal BHF+TBF cal- 
culations we require A = -0.0333 MeV and U = 0.00038 MeV, yielding a 




saturation point at &f ~ 1-36 fm ,E/Am — 15.5 MeV, and an incompress- 
ibility K » 210 MeV. 

These values of A and U have been obtained by using the Argonne vis two- 
body force [20] both in the BHF and in the variational many-body theories. 
However, the required repulsive component (~ U) is much weaker in the BHF 
approach, consistent with the observation that in the variational calculations 
usually heavier nuclei as well as nuclear matter are underbound. Indeed, less 
repulsive TBF became available recently [21] in order to address this problem. 

3. EOS of nuclear matter from different TBF 

In Fig. 3 we compare the different components Vq w , Vji n , V R , Eq. (7), and 
■> V§-> Vc-, Vfi Vt, Eq. (8), of the averaged phenomenological and micro- 
scopic TBF potentials in symmetric matter at normal density. One notes that 
the attractive components V<¥ , V£ n and V™ , Vf roughly correspond to each 
other, whereas the repulsive part (V R vs. Vc) is much larger for the micro- 
scopic TBF. With the choice of parameters A and U given above, one would 
therefore expect a more repulsive behaviour of the microscopic TBF, which is 
indeed confirmed in the following. 

Let us now confront the EOS predicted by the phenomenological TBF and 
the microscopic one. In both cases the BHF approximation has been adopted 
with same two-body force (Argonne v\%). In the left panel of Fig. 4 we display 
the equation of state both for symmetric matter (lower curves) and pure neutron 
matter (upper curves). We show results obtained for several cases, i.e., i) only 
two-body forces are included (dotted lines), ii) TBF implemented within the 
phenomenological Urbana IX model (dashed lines), and iii) TBF treated within 
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Figure 4. Left plot: Binding energy per nucleon of symmetric nuclear matter (lower curves 
using a given linestyle) and pure neutron matter (upper curves), employing different TBF. Right 
plot: Corresponding symmetry energy of nuclear matter. 

the microscopic meson-exchange approach (solid lines). We notice that the 
EOS for symmetric matter with TBF reproduces the correct nuclear matter 
saturation point. Moreover, the incompressibility turns out to be compatible 
with the values extracted from phenomenology, i.e., K m 210 MeV. Up to a 
density of p ps 0.4 fm -3 the microscopic and phenomenological TBF are in 
fair agreement, whereas at higher density the microscopic TBF turn out to be 
more repulsive. 

Within the BHF approach, it has been verified [15, 22] that a parabolic ap- 
proximation for the binding energy of nuclear matter with arbitrary proton frac- 
tion x is well fulfilled, 

|(p,x) « j(p,x = 0.5) + (1 - 2x) 2 E sym (p) , (9) 

where the symmetry energy E sym can be expressed in terms of the difference of 
the energy per particle between pure neutron (x = 0) and symmetric (x = 0.5) 
matter: 

ld(E/A), , E, . E , 
E sym ( P ) = -4^7j^(^°) « -j(P^) ~ ^0.5) . (10) 

In the right panel of Fig. 4 we display the symmetry energy as a function of 
the nucleon density p for different choices of the TBF. We observe results in 
agreement with the characteristics of the EOS shown in the left panel. Namely, 
the stiffest equation of state, i.e., the one calculated with the microscopic TBF, 
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yields larger symmetry energies compared to the ones obtained with the Ur- 
bana phenomenological TBF. Moreover, the symmetry energy calculated (with 
or without TBF) at the saturation point yields a value E syni k, 30 MeV, com- 
patible with nuclear phenomenology. 

4. Neutron star structure 

In order to study the effects of different TBF on neutron star structure, we 
have to calculate the composition and the EOS of cold, catalyzed matter. We 
require that the neutron star contains charge neutral matter consisting of neu- 
trons, protons, and leptons (e~, p~) in beta equilibrium. Using the various 
TBF discussed above, we compute the proton fraction and the EOS for charge 
neutral and beta-stable matter in the following standard way [23, 24]: The 
Brueckner calculation yields the energy density of lepton/baryon matter as a 
function of the different partial densities, 

E 

t{Pn, Pp, Pe, Pn) = (p n m n + P p m p ) + (p n + Pp) — (p n ,p p ) 

1 (37^)5/3 

where we have used ultrarelativistic and nonrelativistic approximations for the 
energy densities of electrons and muons, respectively. The various chemical 
potentials (of the species i = n,p,e,p) can then be computed straightfor- 
wardly, 

de 

IH = 7i- , (12) 

dpi 

and the equations for beta-equilibrium, 

Pi = hpn ~ QiPe , (13) 

(bi and denoting baryon number and charge of species i) and charge neutral- 
ity. 

PiU = . ( 14 ) 

i 

allow to determine the equilibrium composition Pi(p) at given baryon density 
P = pn + Pp and finally the EOS, 

2 d e( Pi (p)) de 

P{P) = P = Pl~~ e = PP™ ~ e ■ ( 15 ) 

dp p dp 




In order to calculate the mass-radius relation, one has then to solve the well- 
known Tolman-Oppenheimer-Volkov equations [23], 

Gm (e + p)(l + Airr 3 p/m) 
r 2 1 — 2Gm/r 

4vrr 2 e , (17) 

with the newly constructed EOS for the charge neutral and beta-stable case as 
input (supplemented by the EOS of Feynman-Metropolis-Teller [25], Baym- 
Pethick-Sutherland [26], and Negele-Vautherin [27] for the outer part of the 
neutron star, p < 0.08 fm~ 3 ). The solutions provide information on the inte- 
rior structure of a star, p(r), as well as the mass-radius relation, M{R). 

The results are shown in Fig. 5. We notice that the EOS calculated with 
the microscopic TBF produces the largest gravitational masses, with the maxi- 
mum mass of the order of 2.3 M©, whereas the phenomenological TBF yields 
a maximum mass of about 1.8 M . In the latter case, neutron stars are char- 
acterized by smaller radii and larger central densities, i.e., the Urbana TBF 
produce more compact stellar objects. For completeness, we also show a se- 
quence of stellar configurations obtained using only two-body forces. In this 
case the maximum mass is slightly above 1.6 M Q , with a radius of 9 km and a 
central density equal to 9 times the saturation value. 

However, these results should be considered as only provisory, since it is 
well known that the inclusion of hyperons [15, 16] or quark matter [28] may 
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strongly affect the structure of the star, in particular reducing substantially the 
maximum mass. We discuss this point now in detail. 

5. Hyperons in nuclear matter 

While at moderate densities p « p the matter inside a neutron star consists 
only of nucleons and leptons, at higher densities several other species of parti- 
cles may appear due to the fast rise of the baryon chemical potentials with den- 
sity. Among these new particles are strange baryons, namely, the A, S, and S 
hyperons. Due to its negative charge, the Yr hyperon is the first strange baryon 
expected to appear with increasing density in the reaction n + n — > p + 
in spite of its substantially larger mass compared to the neutral A hyperon 
(M s - = 1197 MeV, M A = 1116 MeV). Other species might appear in stel- 
lar matter, like A isobars along with pion and kaon condensates. It is therefore 
mandatory to generalize the study of the nuclear EOS with the inclusion of the 
possible hadrons, other than nucleons, which can spontaneously appear in the 
inner part of a neutron star, just because their appearance is able to lower the 
ground state energy of the nuclear matter dense phase. In the following we will 
concentrate on the production of strange baryons and assume that a baryonic 
description of nuclear matter holds up to densities as those encountered in the 
core of neutron stars. 

As we have seen in the previous sections, the nuclear EOS can be calcu- 
lated with good accuracy in the Brueckner two hole-line approximation with 
the continuous choice for the single-particle potential, since the results in this 
scheme are quite close to the full convergent calculations which include also 
the three hole-line contribution. It is then natural to include the hyperon de- 
grees of freedom within the same approximation to calculate the nuclear EOS 
needed to describe the neutron star interior. To this purpose, one requires in 
principle nucleon-hyperon (NY) and hyperon-hyperon (YY) potentials. In our 
work we use the Nijmegen soft-core NY potential [29] that is well adapted to 
the available experimental NY scattering data. Unfortunately, up to date no 
YY scattering data and therefore no reliable YY potentials are available. We 
therefore neglect these interactions in our calculations, which is supposedly 
justified, as long as the hyperonic partial densities remain limited. Also, for 
the following calculations the fig NN potential together with the phenomeno- 
logical TBF introduced previously, are used. 

With the NN and NY potentials, the various G matrices are evaluated by 
solving numerically the Brueckner equation, which can be written in operato- 
rial form as [13, 15] 



G ab [W] = V ab + J2J2 V *' 

C p,p> 



\pp 



Qc 

W -E c + ie 



pp 



G cb [W] 



(18) 
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where the indices a, b, c indicate pairs of baryons and the Pauli operator Q and 
energy E determine the propagation of intermediate baryon pairs. In a given 
nucleon-hyperon channels c = (NY) one has, for example, 



E, 



(NY) 



rriN + my + 



k 2 
N 



+ 



2m 7v 2my 



+ U N (k N ) + U Y (k Y ) 



(19) 



The hyperon single-particle potentials within the continuous choice are given 
by 



U Y (k)=Re E ( kk> G (NY)(NY)[E {N Y)(k,k')} 



and similar expressions of the form 



kk' 



(20) 



u N (k)= u ( N N ' ] (k)+ u N } (k) 



(21) 



N'=n,p 



F=S-,A 



apply to the nucleon single-particle potentials. The nucleons feel therefore 
direct effects of the other nucleons as well as of the hyperons in the environ- 
ment, whereas for the hyperons there are only nucleonic contributions, because 
of the missing hyperon-hyperon potentials. The equations (18-21) define the 
BHF scheme with the continuous choice of the single-particle energies. Due 
to the occurrence of Un and Uy in Eq. (19) they constitute a coupled system 
that has to be solved in a self-consistent manner. In contrast to the standard 
purely nucleonic calculation there is now an additional coupled channel struc- 
ture, which renders a self-consistent calculation quite time-consuming. 

Once the different single-particle potentials are known, the total nonrela- 
tivistic baryonic energy density, e, can be evaluated: 



E 



i=n,p,H ,A 



k F } dk k 2 



k 



1. 



mi + - — + -Ui(k) 



CNN 



Y=£- A 



^dkk 2 



7T^ 



2m 



my + 



2my 



+ U Y a) (k) + U^(k) 



(22) 



(23) 



where e^N is the nucleonic part of the energy density, Eq. (11). Using for 
example an effective mass approximation for the hyperon single-particle po- 
tentials, one could write the last term due to the nucleon-hyperon interaction 

as 



€NY 



E 

Y=£- A 



PY 



my + Uy(0) 



+ 



1 (3vrVy) 5/3 ^ 



2m\ 



5vr 2 



(24) 
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Figure 6. The single-particle potentials of nucleons n, p and hyperons E~, A in baryonic 
matter of fixed nucleonic density pjv = 0.4 fm -3 , proton density p p /pN = 0.2, and vary- 
ing £~ density ps/pN = 0.0, 0.2, 0.5. The vertical lines represent the corresponding Fermi 
momenta of n, p, and £~. For the nucleonic curves, the thick lines represent the complete 
single-particle potentials Un, whereas the thin lines show the values excluding the E~ contri- 



bution, i.e., U ( " ] + U 



Ap) 



which should be added to Eq. (11). 

The different single-particle potentials involved in the previous equations 
are illustrated in Fig. 6, where neutron and proton densities are fixed, given by 
Pn = 0.4 fm -3 and p p /pn = 0.2, and the density is varied. Under these 
conditions the £~ single-particle potential is sizeably repulsive, while Ua is 
still attractive (see also Ref. [15]) and the nucleons are both strongly bound. 
The £~ single-particle potential has a particular shape with an effective mass 
m* /m close to 1, whereas the A effective mass is typically about 0.8 and the 
nucleon effective masses are much smaller. 

The knowledge of the energy density allows then to compute EOS and neu- 
tron star structure as described before, now making allowance for the species 
i = n, p, £~ , A, e~,fi~. The main physical features of the nuclear EOS which 
determine the resulting compositions are essentially the symmetry energy of 
the nucleon part of the EOS and the hyperon single-particle potentials inside 
nuclear matter. Since at low enough density the nucleon matter is quite asym- 
metric, the small percentage of protons feel a deep single-particle potential, 
and therefore it is energetically convenient to create a £~ hyperon, since then 
a neutron can be converted into a proton. The depth of the proton potential is 
mainly determined by the nuclear matter symmetry energy. Furthermore, the 
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Figure 7. The equilibrium composition of asymmetric and /3-stable nuclear matter comprising 
(a) only nucleons, (b) noninteracting hyperons, (c) interacting hyperons. 



potentials felt by the hyperons can shift substantially the threshold density at 
which each hyperon sets in. 

In Fig. 7 we show the chemical composition of the resulting /3-stable and 
asymmetric nuclear matter containing hyperons. We observe rather low hy- 
peron onset densities of about 2-3 times normal nuclear matter density for the 
appearance of the £~ and A hyperons. (Other hyperons do not appear in the 
matter). Moreover, an almost equal percentage of nucleons and hyperons are 
present in the stellar core at high densities. A strong deleptonization of matter 
takes place, since it is energetically convenient to maintain charge neutrality 
through hyperon formation rather than /3-decay. This can have far reaching 
consequences for the onset of kaon condensation. 

The resulting EOS is displayed in Fig. 8. The upper curves show the EOS 
when stellar matter is composed only of nucleons and leptons. The inclusion 
of hyperons (lower curves) produces a much softer EOS, which turns out to 
be very similar to the one obtained without TBF This is quite astonishing be- 
cause, in the pure nucleon case, the repulsive character of TBF at high density 
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Figure 8. The EOS for hyperon-free (solid line) and hyperon-rich (dashed lines) matter. 



increases the stiffness of the EOS, thus changing dramatically the equation of 
state. However, when hyperons are included, the presence of TBF among nu- 
cleons enhances the population of £~ and A because of the increased nucleon 
chemical potentials with respect to the case without TBF, thus decreasing the 
nucleon population. Of course, this scenario could partly change if hyperon- 
hyperon interactions were known or if TBF would be included also for hyper- 
ons, but this is beyond our current knowledge of the strong interaction. 

The consequences for the structure of the neutron stars are illustrated in 
Fig. 9, where we display the resulting neutron star mass-radius curves, com- 
paring now results obtained with different nucleonic TBF, in analogy to Fig. 5. 
One notes that while in Fig. 5 the different TBF still yield quite different max- 
imum masses, the presence of hyperons equalizes the results, leading now to a 
maximum mass of less than 1.3 solar masses for all the nuclear TBF. 

This surprising result is due to the strong softening of the baryonic EOS 
when including hyperons as additional degrees of freedom, and we do not 
expect substantial changes when introducing refinements of the theoretical 
framework, such as hyperon-hyperon potentials, hyperonic TBF, relativistic 
corrections, etc. The only remaining possibility in order to reach larger max- 
imum masses appears the transition to another phase of dense (quark) matter 
inside the star. This will be discussed in the following. 

6. Quark matter 

The results obtained with a purely baryonic EOS call for an estimate of the 
effects due to the hypothetical presence of quark matter in the interior of the 
neutron star. Unfortunately, the current theoretical description of quark matter 
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Figure 9. Neutron star gravitational mass vs. radius (left panel) and the central baryon density 
p c (right panel). Calculations involving different nucleonic TBF are compared. 



is burdened with large uncertainties, seriously limiting the predictive power 
of any theoretical approach at high baryonic density. For the time being we 
can therefore only resort to phenomenological models for the quark matter 
EOS and try to constrain them as well as possible by the few experimental 
information on high density baryonic matter. 

One important condition is constituted by the fact that certainly in symmet- 
ric nuclear matter no phase transition is observed below pa 3po- In fact some 
theoretical interpretation of the heavy ion experiments performed at the CERN 
SPS [30] points to a possible phase transition at a critical density p c pa 6po ~ 
1/fm 3 . We will in the following take this value for granted and use an ex- 
tended MIT bag model [31] (requiring a density dependent bag "constant") 
that is compatible with this condition. 

We first review briefly the description of the bulk properties of uniform 
quark matter, deconfined from the /3-stable hadronic matter mentioned in the 
previous section, by using the MIT bag model [31]. The thermodynamic po- 
tential of / = u, d, s quarks can be expressed as a sum of the kinetic term and 
the one-gluon-exchange term [32, 33] proportional to the QCD fine structure 
constant a s , 
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where m,f and pf are the / current quark mass and chemical potential, respec- 
tively, yf = fif/rrif, Xf = ^Jyj — 1, and <j rcn = 313 MeV is the renormaliza- 
tion point. The number density pf of / quarks is related to Qf via 

and the total energy density for the quark system is written as 

e Q (p u ,p d ,Ps) = J2(tt f + fifPf) + B , (27) 
/ 

where B is the energy density difference between the perturbative vacuum and 
the true vacuum, i.e., the bag "constant." 

In the original MIT bag model the bag constant B pa 55 MeV fm~ 3 is used, 
while values B pa 210 MeV fm~ 3 are estimated from lattice calculations [34]. 
In this sense B can be considered as a free parameter. We found, however, 
that a bag model involving a constant (density independent) bag parameter 
B, combined with our BHF hadronic EOS, will not yield the required phase 
transition in symmetric matter at p c pa 6p pa 1/fm 3 [28]. This can only be 
accomplished by introducing a density dependence of the bag parameter. (The 
dependence on asymmetry is neglected at the current level of investigation). In 
practice we use a gaussian parametrization, 



B(p) = + (B - Boo) exp 



-P 



(28) 



with B^ = 50 MeVfm- 3 , B = 400 MeVfm" 3 , and = 0.17, displayed in 
Fig. 10(a). 

For the description of a pure quark phase inside the neutron star, as for 
neutrino-free baryonic matter, the equilibrium equations for the chemical po- 
tentials, 

Pd = Ps = Pu + Pe , (29) 

must be supplemented with the charge neutrality condition and the total baryon 
number conservation, 

= ^(2p u - p d - p s ) - p e , (30) 

P = ^(.Pu + Pd + Ps) , (31) 

in order to determine the composition pf(p) and the pressure of the quark 
phase, 

P Q (p)=P^-e Q . (32) 
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Figure 10. (a) Bag constant B versus baryon number density, (b) EOS including both 
hadronic and quark components. The shaded region, bordered by two dots, indicates the mixed 
phase (MP) of quarks and hadrons, while HP and QP label the pure hadron and quark phases. 



However, a more realistic model for the phase transition between baryonic 
and quark phase inside the star is the Glendenning construction [35], which 
determines the range of baryon density where both phases coexist. The es- 
sential point of this procedure is that both the hadron and the quark phase are 
allowed to be separately charged, still preserving the total charge neutrality. 
This implies that neutron star matter can be treated as a two-component sys- 
tem, and therefore can be parametrized by two chemical potentials like elec- 
tron and baryon chemical potentials fi e and p, n . The pressure is the same in 
the two phases to ensure mechanical stability, while the chemical potentials of 
the different species are related to each other satisfying chemical and beta sta- 
bility. The Gibbs condition for mechanical and chemical equilibrium at zero 
temperature between both phases reads 

PH(fJ>e, Mn) = Pq{Pc Hn) = Pm(Pu) • (33) 

From the intersection of the two surfaces representing the hadron and the quark 
phase one can calculate the equilibrium chemical potentials of the mixed phase, 
as well as the charge densities pf and p® and therefore the volume fraction x 
occupied by quark matter in the mixed phase, 

P H 

Pc - Pc 

From this, the baryon density pu and the energy density em of the mixed 
phase can be calculated as 



PM = XPQ + (1 - X)PH , 
£M = X£Q + (1 - X)ZH ■ 



(35) 
(36) 
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Figure 11. The gravitational mass (in units of the solar mass M ) versus the normalized 
central energy density (e = 156 MeVfm -3 ) (left panel) and versus the equatorial radius 
(right panel). The thin lines represent static equilibrium configurations, whereas the thick lines 
display configurations rotating at their respective Kepler frequencies. Several different stellar 
matter compositions are considered (see text for details). 

The EOS resulting from this procedure is shown in Fig. 10(b), where the pure 
hadron, mixed, and pure quark matter portions are indicated. The mixed phase 
begins actually at a quite low density around po. Clearly the outcome of 
the mixed phase construction might be substantially changed, if surface and 
Coulomb energies were taken into account [36]. For the time being these are, 
however, unknown and have been neglected. 

The final result for the structure of hybrid neutron stars is shown in Fig. 11, 
displaying mass-radius and mass-central density relations. It is evident that 
the most striking effect of the inclusion of quark matter is the increase of the 
maximum mass, now reaching about 1.5 M & . At the same time, the typical 
neutron star radius is reduced by about 3 km to typically 9 km. Hybrid neutron 
stars are thus more compact than purely hadronic ones and their central energy 
density is larger. For completeness, the figure shows besides static neutron star 
configurations also those rotating at the maximum (Kepler) frequency [37]. In 
that case one observes a further enhancement of the maximum mass to about 
1.8 Mq, and an increase of the typical equatorial radius by about 1 km. 

Finally, in Fig. 12 we display the Kepler periods Pk (= 2tt/Uk) versus the 
rotational star mass for several different stellar sequences based on different 
EOS. Purely hadronic stars, shown by the dotted and long-dashed lines respec- 
tively, show instability against mass shedding first, because of their relatively 
large equatorial radii. Their limiting mass configurations are characterized by 
values of the Kepler period larger than half a millisecond, in agreement with 
results usually found in the literature [38]. In contrast, hybrid stars can reach 
stable periods smaller than half a millisecond. 
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Figure 12. Kepler period versus the rotational mass for purely hadronic stars as well as hybrid 
stars. The following core compositions are considered: i) nucleons and leptons (dotted line); 
ii) nucleons, hyperons, and leptons (dashed line); iii) hadrons, quarks, and leptons (solid line). 
The shaded area represents the current range of observed data. 

7. Conclusions 

In this contribution we reported the theoretical description of nuclear mat- 
ter in the BHF approach and its various refinements, with the application to 
neutron star structure calculation. We pointed out the important role of TBF 
at high density, which is, however, strongly compensated by the inclusion of 
hyperons. The resulting hadronic neutron star configurations have maximum 
masses of only about 1.3 M , and the presence of quark matter inside the star 
is required in order to reach larger values. 

Concerning the quark matter EOS, we found that a density dependent bag 
parameter B(p) is necessary in order to be compatible with the CERN-SPS 
findings on the phase transition from hadronic to quark matter. Joining the 
corresponding EOS with the baryonic one, maximum masses of about 1.6 M G 
are reached, in line with other recent calculations of neutron star properties 
employing various phenomenological RMF nuclear EOS together with either 
effective mass bag model [39] or Nambu-Jona-Lasinio model [40] EOS for 
quark matter. 

The value of the maximum mass of neutron stars obtained according to our 
analysis appears rather robust with respect to the uncertainties of the nuclear 
and the quark matter EOS. Therefore, the experimental observation of a very 
heavy (M > 1.6M ) neutron star, as claimed recently by some groups [41] 
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(M pa 2.2 Mq), if confirmed, would suggest that either serious problems are 
present for the current theoretical modelling of the high-density phase of nu- 
clear matter, or that the assumptions about the phase transition between hadron 
and quark phase are substantially wrong. In both cases, one can expect a well 
defined hint on the high density nuclear matter EOS. 

We would like to thank our collaborators J. Cugnon, A. Lejeune, U. Lom- 
bardo, F. Mathiot, P.K. Sahu, F. Weber, X.R. Zhou, and W. Zuo. 
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